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This paper presents the results of oscillation analysis in K2K experiment. The results show 
indications of neutrino oscillation and give a new constraint on the oscillation parameters. The 
difference of neutrino masses squared Am 2 lies between 1.5 and 3.9xl0 -3 eV 2 at sin 2 20 = 1 
with the confidence level of 90%. In addition to these results, a brief overview of future 
long-baseline neutrino experiment in Japan, JHF-^ experiment, is also given in this paper. 



1 Introduction 



Since a discovery of neutrino oscillation in atmospheric neutrinos by Super-Kamiokande (SK) , 
a lot of attention has been directed to this new phenomena. The zenith angle distribution of 
atmospheric neutrinos observed in SK shows a clear deficit of upward- going Va, which is well 
explained by two- flavor Vp—v T oscillation with Am 2 ~ 3 x 1CP 3 eV 2 and sin 2 29 ~ 1. 

The KEK-to-Kamioka long-baseline neutrino experiment (K2K) uses an accelerator-produced 
neutrinos with a neutrino flight length of 250 km to probe the same Am 2 region as that explored 
by atmospheric neutrinos. Recently, we have performed analysis for neutrino oscillation using 
the information of flux at SK together with the shape of energy spectrum. This analysis is 
based on the data taken from June, 1999 to July 2001 (K2K-I) corresponding to 4.8 x 10 19 pro- 
tons on target (POT), which is approximately half of allocated POT. Since there already exist 
a lot of descriptions on K2K experiment 2 ' 3,4 ' 5 , only the procedure of oscillation analysis and its 
results are presented in this paper. The description in this paper is based on the publication 6 . 

While the first generation neutrino oscillation experiments have successfully being done and 
the oscillations in the atmospheric and solar region are confirmed by several experiments, v^-v e 
oscillation has not been seen except for LSND and is still in mystery. In order to discover the z/«— 



u e oscillation and to open possibilities for detailed studies of neutrino mass ,lepton sector mixing 
and farther possibilities for discovery of new physics, JHF-to-Kamioka long-baseline neutrino 
experiment, JHF-z^, is planned in Japan as a next generation neutrino experiment. A brief 
overview of JHF-z/ experiment is given in the last section. 

2 Strategy for Oscillation Analysis in K2K 

The neutrino beam in K2K experiment is produced by a 12 GeV proton beam from the KEK 
proton synchrotron. After protons hit an aluminum target, the produced positively charged 
particle, mainly 7r + , are focused by a pair of pulsed magnetic horns 7 . The pion momentum and 
angular distributions downstream of second horn are occasionally measured by gas-Cherenkov 
detector (PIMON) 8 in order to verify the beam Monte Carlo (MC) simulation and to estimate 
the errors on the flux prediction at SK 2 . The neutrino beam produced from the decays of these 
particles in a decay pipe of 200 m long are 98% pure muon neutrinos with a mean energy of 1.3 
GeV. Near neutrino detectors (ND) are located at 300 m from the proton target. ND measure 
the stability and direction of the neutrino beam as well as neutrino flux and energy spectrum 
just after the production. SK is located at 250 km from KEK. Here the neutrino flux and energy 
spectrum after neutrinos travel the length of 250 km are measured. The measurements at ND 
are extrapolated to SK by multiplying far-to-near flux ratio (F/N). The detailed description 
for F/N can be found in the reference 2 . 

In the case of the presence of neutrino oscillation, both the reduction in the number of 
neutrino events and the distortion in the energy spectrum are expected after neutrinos travel a 
fixed flight length, because the oscillation probability depends on neutrino energy: 

P(^-^)=sin 2 20sin 2 ^^, (1) 

where 6 is mixing angle, Am 2 is the difference of neutrino mass squared, E v is neutrino energy, 
and L is the neutrino flight length which is fixed to be 250 km in the K2K case. Therefore both 
the number of observed events and the spectral shape information at SK are used to compare 
them with the expectation from ND measurements. In this analysis, all the beam induced neu- 
trino events observed in the fiducial volume of SK are used to measure the overall suppression 
of flux, while in order to enhance the fraction of charged-current (CC) quasi-elastic (QE) inter- 
actions (y^ + n — > /x + p), 1-ring //-like events (lR/i) are used to study the spectral distortion. 
Since the proton momentum in the CCQE events is typically below Cherenkov threshold, only 
muon is visible. Assuming QE interaction in lR/i and neglecting Fermi momentum, the neutrino 
energy can be reconstructed as 

^nEu + m 2 
m N - E^+p^ cos Ofj, ' 

where tun, E^, m M , p^ and 0^ are nucleon mass, muon energy, muon mass, muon momentum 
and muon scattering angle with respect to the direction of neutrino beam, respectively. 

3 Measurements of Neutrino Flux and Spectrum at ND 

The ND consist of two detector systems: a 1 kilo-ton water Cherenkov detector (1KT) and 
a fine-grained detector (FGD). The flux normalization is measured by 1KT to estimate the 
expected number of events at SK. Since 1KT has the same detector technology as SK, the most 
of the systematic uncertainties on the measurements, mainly on the detection efficiency and 
cross-section of interactions, are canceled in comparison between ND and SK measurements. 



(2) 



SciFi 2 track cos(A@ p ) distribution 

600 
500 



400 
300 
200 
100 



"-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

COS(A0 p ) 

Figure 1: The distribution of cosA# p for 2-track sample in SciFi. The schematic explanation for A9 P is also 
drawn. The events with A0 P < 25 deg. are selected as QE enhanced sample, while those with A9 P > 30 deg. are 

selected as non-QE enhanced sample. 



The energy spectrum is measured by analyzing the muon momentum and angular distributions 
in the both detector systems. The 1KT has high efficiency to reconstruct muons with the 
momentum of below 1 GeV/c and full 4ir coverage in solid angle. However, the 1KT has little 
efficiency to reconstruct muons with the momentum of above 1.5 GeV/c since they exit the 
detector. On the other hand, the FGD has high efficiency to measure the muons above 1 
GeV/c, the relevant energy range can be covered by these two complementary detectors. 



3.1 1KT Analysis 

In the 1KT analysis, the interactions are measured by detecting Cherenkov light emitted 
from produced charged particles. The vertices, direction and momentum of Cherenkov rings are 
reconstructed with the same method as in SK 9 . The event selection criteria for the measurement 
of integrated flux are the same as those in the reference 2 : (a) Deposited energy is greater than 
100 MeV. (b) The reconstructed vertex is inside the 25 t fiducial volume which is defined by a 
2 m radius, 2 m long cylinder oriented to the beam axis, in the upstream side of the detector. 
The expected number of events with the vertices fully contained (FC) in SK fiducial volume is 
estimated to be 8O.1I54. The correlation between energy bins from the spectrum measurement 
at ND and F/N are taken into account in the estimation of the errors, described in detail below. 
The major contribution to the errors come from the uncertainties in the F/N (^'q%) and the 
normalization (±5%). The normalization error is dominated by uncertainties of the fiducial 
volumes due to the vertex reconstruction both in 1KT (4%) and in SK (3%). For the spectrum 
measurement, we impose further cuts in order to select lUfi events: (c) The event has one ring, 
(d) The particle stops inside the 1KT. Among the events after the cuts (a) and (b), 53% have 
one ring. For criteria (d), we require the maximum charge of hit photomultiplier tubes (PMT) 
to be less than 200 p.e. This eliminates effectively the events with muon exiting the 1KT; 68% 
of one ring events retain after this cut. Using this sample, we make the muon momentum and 
angular distributions to estimate the spectrum shape of neutrino, which is described in detail 
below. The largest systematic uncertainty for the spectrum measurement in the 1KT is that on 
energy scale. We estimate the error on energy scale using cosmic-ray muons and beam-induced 
7T° sample and quote ^3% for this error. 




Table 1: The central values of the flux re- weighting parameters for the spectrum fit at ND ($nd) and the 
percentage size of the energy dependent systematic errors on $nd, F/N, and esk- The re weighting parameters 

are given relative to the 1.0-1.5 GeV energy bin. 
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3.2 FGD Analysis 

The FGD consist of a scintillating fiber tracker with water target (SciFi) 3 , plastic scintillator 
trigger counters upstream and downstream of SciFi (TRIG), a lead-glass calorimeter (LG), and 
a muon range detector (MRD) 4 . In the FGD analysis, we use the events which have one or two 
tracks with the vertex within the 5.9 t fiducial volume of SciFi, which is defined as a rectangle 
of 2.2 m x 2.2 m in x and y, and 1 st to 17 th water containers in the beam direction (z). The 
track finding efficiency is 70% for a track passing through three layers of SciFi and close to 100% 
for more than five layers 10 . Three layers is the minimum track length required in this analysis. 
Events which have at least one track passing into MRD are chosen in order to select ^-induced 
CC interactions. Furthermore since SciFi and MRD do not have fine timing information to select 
on-spill events, it is required that the hit in downstream TRIG which matches with SciFi track 
should be in spill timing. The momentum of each track is measured by its range through the 
SciFi, TRIG, LG and MRD with the accuracy of 2.7%. The second track also be reconstructed if 
the proton produced in the QE interaction has a momentum of typically greater than 600 MeV/c. 
In the case where the second track is visible, the kinematical information is used to enhance the 
fraction of QE and non-QE events in the sample. The direction of second track can be predicted 
from the momentum of the first track by assuming QE interaction. The distribution of cosine 
of the angular difference between the predicted and observed second track (cosA# p ) is shown 
in Fig. 1. QE enhanced and non-QE enhanced samples are selected by requiring A9 P < 25 deg. 
and A9 p > 30 deg., respectively. The fraction of the QE events in the QE sample is estimated 
to be 62% by MC simulation, while 82% of events in non-QE sample is estimated to come from 
interactions other than QE. The SciFi events are divided into three event categories; 1-track, 
2-track QE enhanced, and 2-track non-QE enhanced samples. We make muon momentum and 
angular distributions for each event category to estimate the neutrino energy spectrum at ND, 
which is described below. 

3.3 Neutrino Spectrum at ND 

In order to estimate the neutrino energy spectrum at ND, i.e. just after the production, the 
2-dimensional distributions of muon momentum and angle with the respect to the neutrino beam 
direction for four event categories are used: lR/i sample in 1KT, 1-track, 2-track QE enhanced, 
and 2-track non-QE enhanced samples in FGD. A x 2 -fitting method is used to compare these 
data against the MC expectation. The neutrino spectrum is divided into eight energy bins as 
defined in Table 1. The flux in each energy bin is re- weighted relative to the value of the beam 
MC simulation. These re-weighting parameters are normalized such as the bin of E v = 1.0- 



1kt Muon Momentum (FV25t 1-ring u-like FC) 1 kt Muon Angle (FV25t 1-ring u-iike FC) 




Figure 2: The momentum and angular distributions for 1KT lKfi sample. The left figure is the distribution of 
momentum and the right one is that of angle with respect to the beam direction. The hatched histograms are 
the contribution from CCQE interaction in MC simulation. 
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Figure 3: The momentum and angular distributions for SciFi event samples. The left figures are the distributions 
of momentum and the right ones are those of angle with respect to the beam direction for 1-track, 2-track QE 
enhanced, and 2-track non-QE enhanced samples in order of the top, middle, and bottom figures. Here again, 
the hatched histograms are contribution from CCQE interaction in MC simulation. 
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Figure 4: The neutrino energy spectrum at near site inferred from ND data along with the result of beam MC 
simulation. Points with error bars are the resulting spectrum of ND fitting, while histogram is the result of beam 

MC simulation. 



1.5 GeV to be 1.0. An overall normalization is introduced as a free parameter in the fit. The 
parameter -R nqe is used to re- weight the ratio between the QE and non-QE cross-section relative 
to the MC simulation. The systematic uncertainties in the ND measurements, such as the energy 
scales, the track finding efficiencies, the detector thresholds, and so on, are incorporated into 
the fitting. In addition, PIMON is used in order to constrain the re-weighting parameters for 
E v >1.0 GeV. 

The fitting is successfully done, and the x 2 value at the best-fit point is 227.2/197d.o.f. All 
the parameters including the detector systematics are found to lie within their expected errors. 
The best-fit values of the flux re- weighting factors are shown in Table 1. The muon momentum 
and angular distributions for lR/i sample in 1KT and for 1-track, 2-track QE enhanced and 
2-track non-QE enhanced samples are shown in Fig. 2 and 3 with the re-weighted MC overlaid. 
The fitting results are in good agreement with the data. The resulting neutrino spectrum 
inferred from the ND data is shown in Fig. 4 along with the beam MC simulation result. As 
can be seen, the results of beam MC excellently agrees with ND measurements. The errors on 
this measurement are provided in the form of an error matrix, and the correlations between the 
parameters are taken into account in the oscillation matrix using this matrix, later described in 
detail. The diagonal elements of the error matrix are shown in Table 1. 

3.4 Neutrino Interaction Models 

The uncertainty due to neutrino interaction models is studied separately. In QE interaction, 
the axial vector mass in the dipole formula is set to be 1.1 GeV/c 2 as a central value and is 
varied by ±10%. In single pion production, the central value of the axial vector mass is set 
to be 1.2 GeV/c 2 and varied by ±20% 11 . For coherent pion production, two different models 
are compared: One is the Rein and Sehgal model 12 and the other is a model by Marteau et 
al. 13 . For deep inelastic scattering, GRV94 14 and the corrected structure function modeled by 
Bodek and Yang 15 are studied. Marteau model for coherent pion production and Bodek and 
Yang structure function for deep inelastic scattering are employed in this analysis. The choice 
of models causes the fitted value of -R nqe to change by ~20%. Therefore, the error of 20% on 
Rnqe is added to the error matrix. It is found that the choice of models and axial vector mass 
does not affect on the result of spectrum measurement beyond the size of the fitted errors. 
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Figure 5: The distribution of AT(= Tsk — Trek — Ttof) f° r FC events in SK. The upper figure shows the 
distribution of ±500 (isec time window for each reduction step in SK. The lower figure show the distribution of 

±5 /xsec time window for FCFV events. 



4 Far-to-near Spectrum Ratio 

The F/N ratio from the beam MC simulation is used to extrapolate the measurement at ND 
to those at SK. The detailed description about F/N can be found in the reference 2 . The errors 
on F/N are also provided in the form of an error matrix. The errors including the correlations 
on F/N for E u > 1 GeV are estimated based on the PIMON measurements, where PIMON 
is sensitive in this energy region. Since PIMON is insensitive in the region of E u < 1 GeV, 
the errors and correlations are estimated based on the uncertainties in the hadron production 
models used in K2K beam MC simulation. The errors on F/N of E u > 1 GeV and on that of 
E u < 1 GeV are treated such that these are not correlated. The diagonal elements in the error 
matrix are summarized in Table 1. 



5 Event Selection in SK 

The selection criteria for SK events are all same as those described in the reference 2 . Accelerator- 
produced neutrino interaction at SK are selected by comparing two Universal Time Coordinated 
(UTC) time stamps from the global positioning system (GPS), Trek for the KEK-PS beam spill 
start time and Tsk for the SK trigger time. The time difference between two UTC time stamps, 
AT = Tsk — Tkek — Ttof, where Ttof is the time of flight of neutrinos from KEK to SK, 
should be distributed around the interval from to 1.1 /usee to match the width of the beam 
spill of KEK-PS. The condition —0.2 < AT < 1.3 /isec is imposed as a timing cut. Other 
criteria are: (a) There is no detector activity within 30 //sec before the event, (b) The total 
collected p.e. in a 300 nsec time window is greater than 200 (~ 20 MeV deposited energy), (c) 
The number of PMTs in the largest hit cluster in the outer-detector is less than 10. (d) The 
deposited energy is greater than 30 MeV. Finally, a fiducial cut is applied to select only the 
events with fitted vertices inside 22.5 kt fiducial volume, which volume is same as that used in 
SK atmospheric neutrino analysis. The detection efficiency of this selection is 93% for charged- 



current interactions and 68% for neutral-current inelastic interactions, for a total of 79%. Fig. 5 
shows the AT distribution at the various stages of the reduction. A clear peak in time with the 
neutrino beam from KEK-PS is observed in the analysis time window. After the above selection, 
56 FC events are observed in the fiducial volume (FCFV events), while the number of FCFV 
events expected by 1KT measurement is 8O.II5' 2 as described in Sec. 3.1. With the timing cut, 
the expected number of atmospheric neutrino background is approximately 10~ 3 event, which 
is negligible. 

6 Oscillation Analysis 

A two-flavor oscillation analysis, with i/„ disappearance, is performed with the use of the 
maximum-likelihood method. In this analysis, both the number of FCFV events and the energy 
spectrum shape for lR/x are used. The likelihood is defined as 

Aot = •'—norm (N ohs , N cxp (Am 2 , sin 2 26, /)) x £ shapc (Am 2 , sin 2 26, f) x £ sys (/). (3) 

The normalization term £ nor m is the Poisson probability to observe iV bs events when the ex- 
pected number of events is N exp (Am 2 , sin 2 26, /). The spectrum shape term £ sna pe is the product 
of the probability for each 1R;U event to be observed at E£ ec = Ef. 

AW 

Anape = II P ( E ^ Am 2 , sin 2 26, /), (4) 
i=i 

where P is the normalized ££ ec distribution estimated by MC simulation and Ni^ is the number 
of IR/j, events. The symbol / denotes a set of parameters which is constrained by the systematic 
errors. These parameters consist of the re- weighted neutrino spectrum measured by ND ($nd) 5 
the F/N ratio, the reconstruction efficiency of SK (esx) for lR/x events, the re-weighting factor 
for the QE/non-QE ratio R Qqe , the energy scale of SK, and the overall normalization. The errors 
on first three items are energy dependent and correlated between each energy bin. The diagonal 
elements of their error matrices are summarized in Table 1. The error on SK energy scale is 
3%. The parameters / is treated as fitting parameters with an additional constraint term £ sys 
in likelihood". 

In the oscillation analysis, the whole data sample is used in the normalization term £ norm , 
i.e. Nobs = 56. In the spectrum shape term £ s hape, the data taken in June 1999 are discarded 
since the target radius and horn current, hence the spectrum shape in June 1999 were different 
from the rest of running period. The discarded data correspond to 6.5% of total POT. The 
number of 1R/U events used in this analysis is 29, while the number of 1R^ events from the MC 
simulation in the case of no oscillation is 44. 

The likelihood is calculated at each point in the Am 2 and sin 2 26 space to search for the 
point where the likelihood is maximized. The best-fit point in the physical region of oscillation 
parameter space is found to be at (sin 2 26, Am 2 ) = (1.0,2.8 x 10~ 3 eV 2 ). If the whole space 
including unphysical region is allowed, these values are found to be (sin 2 26, Am 2 ) = (1.03, 2.8 x 
10~ 3 eV 2 ). At the best-fit point in the physical region, the total number of predicted events 
is 54.2, which is in agreement with the observation of 56 events within statistical error. The 
observed E r u cc distribution for lR/x sample is shown in Fig. 6 together with the expectation for 
the best-fit oscillation parameters and that for no oscillation. The Kolgomorov-Smirnov (KS) 
test is done to check the consistency between the observed and best-fit ££ ec spectrum. A KS 
probability is obtained to be 79%, and the best-fit spectrum shape agrees with the observation. 



"Another approach to treat the systematic parameters / is performed, which also gives consistent results with 
those described in the main text 6 . 
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Figure 6: The reconstructed E v distribution for the lR/i sample (from method 1). Points with error bars are 
data, while the box histogram is the expected spectrum in the case of no oscillation, where the height of the box 
is the systematic error. The solid line is the best-fit spectrum in the case where the oscillation occurred. These 
histograms are normalized by the number of observed lR/x events (29). In addition, the dashed line shows the 
expectation with no oscillation normalized to the expected number of lR/i events (44). 

The likelihood ratio of the no oscillation case to the best-fit point is computed to estimate 
the probability that the observations are only due to a statistical fluctuation instead of neutrino 
oscillation. The no oscillation probability is calculated to be 0.7%. When only normalization 
or shape information is used, the probabilities are 1.3% and 16%, respectively. The oscillation 
parameters preferred by the total flux suppression and the energy distortion alone also agree 
well with each other (Fig. 7 (right)). The allowed regions of oscillation parameters are evaluated 
by calculating the likelihood ratio of each point to the best-fit point. 

Finally the uncertainties of neutrino interaction models are studied in the same way as 
the spectrum measurement at ND. It is found that the effects of difference between interaction 
models are negligible for all the results since the cancellation is caused by using the same models 
in ND and SK. 

The contours for confidence level of 68%, 90%, and 99% are drawn in Fig. 7 (left). The behav- 
iors of — A(ln£tot) m the case where normalization-only, shape-only, and normalization+shape 
information is used in the analysis are plotted in Fig. 7 (right). This shows that the reduction in 
the neutrino flux and the distortion in energy spectrum abserved in K2K experiment prefer the 
same Am 2 region with each other, as described above. — A(ln£ to t) as a function of sin 2 29 at 
the best-fit Am 2 (= 2.8 x 10~ 3 eV 2 ) and as a function of Am 2 at the best-fit sin 2 29 (=1.0) are 
also shown in Fig. 8. The 90% C.L. region of Am 2 is 1.5-3.9xl0~ 3 eV 2 at sin 2 29 = 1.0. The 
allowed region of oscillation parameters obtained by Super-Kamiokande atmospheric neutrino 
observations are also drawn in Fig. 7. The K2K and Super-Kamiokande results are in good 
agreement. 

7 The Upgrade of K2K Near Detector — SciBar Detector — 

In the case where the oscillation parameter is Am 2 ~ 3 x 10~ 3 eV 2 , the neutrino energy for 
oscillation maximum is ~ 0.6 GeV under the condition of fixed flight length of 250 km. Therefore, 
precise measurement of neutrino spectrum and the knowledge of neutrino interaction in low 
energy region are indispensable in order to maximize the K2K sensitivity. However, low energy 
particles of which momentum is below Cherenkov threshold are invisible in the 1KT, while SciFi 
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Figure 7: Left: Allowed region of oscillation parameters with the confidence level of 68%, 90%, and 99%. Solid 
lines are the results from the K2K oscillation analysis. In addition, the results of the observation of atmospheric 
neutrinos by Super-Kamiokande (dashed lines) are overlaid. Right: The Behaviors of — A(ln£tot) as a function 
of Am 2 along sin 2 26 = 1.0 axis in the case where normalization-only, shape-only, and normalization+shape 
information is used in the analysis. The Am 2 regions indicated by normalization-only and shape-only agrees well 

with each other, region as analysis 




Figure 8: The behaviors of — A(ln£ to t) as a function of sin 26» at Am 2 = 2.8 x 10" 3 eV 2 (left) and as a function 

of Am 2 at sin 2 26 = 1.0 (right). 
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Figure 9: The schematic view of SciBar detector 

have small efficiency for muons below 0.8 GeV/c. 

As an upgrade of K2K near detector, SciBar detector is designed/developed to detect such 
low energy particles with high efficiency. The schematic view of SciBar detector is drawn in 
Fig. 9. The detector consists of extruded plastic scintillators which dimensions are 1.3 cm thick, 
2.5 cm wide, and 3 m long. Scintillators are aligned in X or Y to construct a layer (X-layer 
or Y-layer). A pair of X- and Y-layer is called as a module, and 64 modules are aligned in 
z-direction, or beam direction. A hole of diameter of 1.8 mm is opened at the center of cross 
section along through each scintillator. A wavelength shifting fiber is put in this hole to readout 
the scintillation light from each scintillator. The light transmitted outside the scintillator is read 
by 64ch multi-anode photomultiplier tubes at one end of fibers. 

The main feature of SciBar detector is that SciBar is a full-active tracker in which the 
scintillators play roles of both neutrino interaction target and particle detector at the same 
time. Therefore, SciBar has high efficiency for short tracks, can detect low energy protons down 
to 350 MeV/c. SciBar also measure dE/dx of each particle. Using this information, particle 
identifications and momentum reconstruction can be performed. 

The following study can be done with the use of this excellent detector: Precise measurement 
of neutrino spectrum in low energy. Detailed study of neutrino interactions. Especially non-QE 
interactions are the background against neutrino energy reconstruction in CCQE, therefore to 
understand these background is very important to maximize the oscillation sensitivity. 

Almost all the R&D items have been finished and the basic performances of the detector 
have been studied. Installation work will be done during the summer in 2003 and will be finished 
by the end of September 2003. The data taking will be started from the beam time of October 
2003. 

8 Overview of JHF-f Experiment 

The JHF to Kamioka neutrino project is a second generation long baseline neutrino oscillation 
experiment (JHF-f experiment) that probes physics beyond the Standard Model by high preci- 




Figure 10: The baseline of JHF to Kamioka neutrino project (left) and the layout of JHF (right). 



sion measurements of the neutrino masses and mixing. A high intensity narrow band neutrino 
beam is produced by secondary pions created by a high intensity proton synchrotron at J-PARC 
(JAERI). The neutrino energy is tuned to the oscillation maximum of ~1 GeV for a baseline 
length of 295 km toward the world largest water Cherenkov detector, Super-Kamiokande (left 
figure of Fig. 10). Its excellent energy resolution and particle identification enable us to recon- 
struct the initial neutrino energy, which is compared with the narrow band neutrino energy, 
through the quasi-elastic interaction. 

The JHF-f experiment consists of two phases, phase-I and phase-II. The physics goal of the 
phase-I is a search for the v„ — > v e appearance with a factor of 20 higher sensitivity (sin 2 29^ ~ 
0.5 sin 2 2^i3 > 0.003), the — > v T oscillation measurement with an order of magnitude better 
precision (^(Am 2 ^) = 10~ 4 eV 2 and 5(sin 2 2#23) = 0.01) than existing measurements, and a 
confirmation of the ia, — ► v T oscillation or discovery of sterile neutrinos (v s ) by detecting the 
neutral current events. In the phase-II, an upgrade of the accelerator from 0.75 MW to 4 MW 
in beam power and the construction of 1 Mt Hyper-Kamiokande detector at Kamioka site are 
envisaged. By these upgrade, it is expected that improvement of another order of magnitude in 
the Uu — > v e oscillation sensitivity, a sensitive search for the CP violation in the lepton sector 
(CP phase 5 down to 10°— 20°), and improvement of an order of magnitude in the sensitivity of 
the nucleon decay search. 

The descriptions in this paper are all based on the "Letter of Intent for JHF Neutrino 
Experiment" 16 . Here only the items related to neutrino beam of JHF-z^ experiment and the 
search for v e appearance. More detailed descriptions are found in the Letter of Intent. 

8.1 Neutrino Beam at JHF 

The layout of JHF is drawn in Fig. 10(right). The proton beam is fast-extracted from the 
50 GeV PS in a single turn and transported to the production target. The design intensity of 
the PS is 3.3 x 10 14 protons per pulse (ppp) at a repetition rate of 0.285 Hz (3.5 sec period). The 
resulting beam power is 0.75 MW (2.64 MJ/pulse). The spill width is ~ 5.2 //sec. We define a 
typical one year operation as 10 21 POT, which is corresponding to about 130 days of operation. 
The protons are extracted toward inside of the PS ring, and are bent by 90° to SK direction in 
the transport line using superconducting magnets. The secondary pions (and kaons) from the 
target are focused by electromagnetic horns 17 , and decay in the decay pipe. The length of the 
decay pipe from the target position is 130 m. The first near detectors are located at 280 m from 
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Figure 11: (a) Neutrino energy spectra of charged current interactions. Thick solid, dashed and dash-dotted 
histograms are OAl°, OA2° and OA3°, respectively, (b) Comparison of v e and spectra OA2°. Solid (black) 
histogram is and dashed (red) one is v e . Hatched area is contribution from K decay. The low energy v e 

component is due to fi decay. 



the target. There also exists a plan to put second near detectors (intermediate detectors) at 
2 km from the target. 

The beam configuration is off-axis beam (OAB) to produce a narrow neutrino energy spec- 
trum 18 . The optics is almost same as the usual wide band beam, which is also used in K2K 
experiment, however, the axis of the beam optics is displaced by a few degrees from the direc- 
tion of far detector (off-axis). A finite decay angle selected, the neutrino energy becomes almost 
independent of parent pion momentum as a consequence of the characteristics of the Lorenz 
boost, which provides the narrow spectrum. The peak energy of neutrino can be adjusted by 
choosing the off-axis angle. The distance between JHF and SK is 295 km. The neutrino energy 
will be tuned to between 0.4 and 1.0 GeV, which corresponds to Am 2 between 1.6-4xl0 -3 eV 2 
suggested by recent SK and K2K results 1,6 . Fig. 11(a) shows expected neutrino energy spectrum 
of charged current interactions at SK in MC simulation. The OAB is roughly a factor of 3 more 
intense than usual narrow band beam which is made by selecting the momentum of pions using 
dipole magnet placed between two horns. Fluxes and numbers of interactions are summarized in 
Table 2. 

The v e contamination in the beam is expected to be 1% at the off-axis angle of 2° (OA2°). 
The sources of v e are ir — > fj, — > e decay chain and K decay {K e ^). Their fractions are fi- 
decay: 37%, K-decay: 63% for OA2°. The energy spectra of the u e contamination are plotted 
in Fig. 11(b). At the peak energy of the i/„ spectrum, the fe/f^ ratio is as small as 0.2% in 
OAB. This indicates that beam u e background is greatly suppressed (factor ~ 4) by applying 
an energy cut on the reconstructed neutrino energy. 

8. 2 Search for v e Appearance 

The JHF neutrino beam has small v e contamination (0.2% at the peak energy of OAB) and the 
v e appearance signal is enhanced by tuning the neutrino energy to be its expected oscillation 
maximum. Therefore, JHF-f experiment has an excellent opportunity to discover v e appearance 
and hence to measure #13. The sensitivity on v e appearance is studied based on the full MC 
simulations and analysis of SK and K2K experiments. 

The v e appearance signal is searched for in the CCQE interaction, for which the energy 



Table 2: Summary of beam simulation. The peak energy -E pea k is in GeV. The flux is given in 

10 6 /cm 2 /yr, 

and 

the Vejv^ flux ratio is in %. The ratio in the "total" column is the one integrated over neutrino energy and the 
column "iSpcak" is the ratio at the peak energy of spectrum. The normalization for the number of interactions 
are /22.5kt/yr. The numbers outside (inside) the bracket are number of total (CC) interactions. 



Beam 


Speak 


Flux 


Ve/"n (%) 
total Speak 


# of interactions 


OA2° 
OA3° 


0.7 
0.55 


19.2 0.19 
10.6 0.13 


1.00 0.21 
1.21 0.20 


3100(2200) 60(45) 
1100( 800) 29(22) 



Table 3: Number of events and reduction efficiency of "standard" 1-ring e-like cut and n cut for 5 year exposure 
(5 x 10 21 POT) OA2°. For the calculation of oscillated v e , Am 2 = 3 x 10~ 3 eV 2 and sin 2 26> Me = 0.05 is assumed. 



OAB 2° 


Vn CC NC Beam v e Oscillated v e 


(1) Generated in F.V. 

(2) 1R e-like 

(3) e/ir° separation 

(4) 0.4 GeV < E T U CC < 1.2 GeV 


10713.6 4080.3 292.1 301.6 
14.3 247.1 68.4 203.7 
3.5 23.0 21.9 152.2 
1.8 9.3 11.1 123.2 



of neutrino can be calculated using kinematics. Since the proton momentum from the QE 
interaction is usually below the Cherenkov threshold, the signal has only a single electro-magnetic 
shower ring (IRe). 

The standard SK atmospheric neutrino analysis criteria are used to select IRe events: single 
ring, electron like (showering), visible energy greater than 100 MeV, and no decay electrons. 
Reduction of number of events by the "standard" 1-ring e-like cut for charged and neutral 
current events are listed in Table 3. The excellent e/fi separation capability and fi — > e de- 
tecting capability are key features of the effective elimination of charged current and all of 
the inelastic events which contain charged it. The remaining background events at this stage 
are predominantly from single 7T° production through neutral current interactions and from v e 
contamination in the beam. In order to reduce the background events which come from 7r°, 
further cuts are applied based on following information; (1) angle between v beam direction 
and e-like ring, (2) invariant mass of two photons assuming that the event contains two rings 
which cannot be separated by "standard" algorithm, (3) difference between single and double 
ring likelihood, (4) energy fraction of lower energy ring, E(jf)+E(j 2 ) ' assumrn § the event con- 
tains two rings. Table 3 lists the number of events after this e/ir separation. Extra rejection 
of an order of magnitude (23/247.1) in the v„ neutral current background is achieved with the 
signal acceptance of 152.2/203.7=75%. Fig. 12 (left) shows the reconstructed neutrino energy 
distributions for 5 years. The oscillation parameters of Am 2 = 3 x 10~ 3 eV 2 and sin 2 2#i3 = 0.1 
are assumed. A clear appearance peak is seen at the oscillation maximum of E v ~ 0.75 GeV. 
The right plot of Fig. 12 show 90% and 3a limits as a function of the years of operation with 
the systematic uncertainty of background subtraction to be 2%, 5% and 10%. The sensitivity 
of sin 2 2#i3 = 0.006 at 90% confidence level can be achieved in five years of operation. Fig. 13 
shows 90%C.L. contours of sensitivity for 5 year exposure assuming 10% systematic uncertainty 
in background subtraction. 

9 Summary 

In the K2K experiment, both the number of observed neutrino events and the observed energy 
spectrum at SK are consistent with neutrino oscillation. It is concluded that the probability 
that our measurements are explained by statistical fluctuation without oscillation is less than 




Figure 12: (Left) Expected reconstructed neutrino energy distributions of expected signal+BG, total BG, and 
BG from interactions for 5 years exposure of OA2°. (Right) Expected 90%CL sensitivity (thick lines) and 
3<r discovery contours (thin lines) as the functions of exposure time of OA2°. In left figure, expected oscillation 
signals are calculated with the oscillation parameters of Am 2 = 3 x 10 -3 eV 2 and sin 2 26* Me (a effective mixing 
angle = sin 2 #23 • sin 2 2^13) = 0.05. In the right figures, three different contours correspond to 10%, 5%, and 2% 

uncertainty in the background estimation. 
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Figure 13: The 90% C.L. sensitivity contours in 5 years of operation with OA2°. The 90% C.L. excluded region 
of CHOOZ is plotted as a comparison. sm 2 #23 is assumed to be 0.5 and the possible contribution due to #12 term 
is assumed to be small compared to the one due to #13 term. 



1%. The best-fit point in the oscillation parameter space is (sin 2 29, Am 2 ) = (1.0,2.8 eV 2 ) 
and the 90% C.L. region of Am 2 at sin 2 26 = 1.0 is 1.5-3.9 x 10~ 3 eV 2 , which agrees well 
with the atmospheric neutrino observations by Super-Kamiokande. After the accident of Super- 
Kamiokande in November 2001, the reconstruction work have been done by the beginning of 
December 2002, and K2K is now running, as K2K-II, to collect 2 times more data, and finally 
up to planned ~ 10 20 POT. As the upgrade of K2K near detectors, SciBar detector will fully 
installed in summer of 2003. The SciBar detector will enable us to study lower energy neutrino 
interactions and cross-section of neutrino interactions. 

In order to confirm precisely the results of the first generation neutrino experiments and to 
explore the physics beyond the Standard Model, JHF to Kamioka neutrino experiment, JHF-f 
experiment, is planned. With the use of high intensity neutrino beam, z/„ — > v e oscillation will 
be explored down to sin 2 26 13 = 0.006 at Am 2 ~ 3 x 10" 3 eV 2 with 90% C.L., which is 20 times 
higher sensitivity than CHOOZ experiment. The R&D has already started and the experiment 
is planned to start in 2007. 
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